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Abstract 

•rH , 

^ . We consider the lepton sector of the Standard Model and allow for an arbitrary 

number of Higgs doublets and, moreover, for the presence of right-handed neutrino 
singlets which enable the seesaw mechanism. In this framework, we identify and 
calculate the dominant one-loop radiative corrections to the tree-level mass ma- 
trix of the light neutrinos. The interesting feature is that both the tree-level and 
the one-loop contributions to the light-neutrino mass matrix are quadratic in the 
Yukawa couplings, with the effect that the one- loop contribution is smaller than the 
tree-level one mainly because of the one-loop factor (167T 2 ) . We also point out 
the possibility of generating radiatively — in this framework — the ratio of solar over 
atmospheric neutrino mass-squared differences, as needed for the large-mixing-angle 
MSW solution of the solar-neutrino problem. 
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1 Introduction 



It is believed that the final confirmation of small but non-zero neutrino masses is around 
the corner, since neutrino oscillations |J provide an excellent and natural solution for 
the atmospheric- and solar-neutrino problems. In the latter case, matter effects in the 
neutrino oscillations seem to play a crucial role. For reviews, see for instance Ref. ||; 
for the latest solar-neutrino results, see Ref. Q and the references therein. 

The results of the solar- and atmospheric-neutrino experiments have brought about 
an upsurge of model building for the lepton masses and mixings — for reviews of lepton- 
mass-matrix textures and models see Ref. ||. In this effort, the seesaw mechanism H] 
figures prominently as a means for obtaining small Majorana masses for the neutrinos. 

In this paper we discuss a very simple extension of the Standard Model (SM): we add 
to the SM right-handed neutrino singlets, and moreover admit an arbitrary number of 
Higgs doublets. The first addition allows to incorporate the seesaw mechanism in the 
SM. As shown in Ref. 0, this simple extension of the SM provides a framework in which 
interesting models can be constructed, which reproduce in a natural way at least part 
of the striking features of the neutrino masses and mixings. This extension of the SM 
has two large mass scales: the electroweak scale, of about 100 GeV, which is the order 
of magnitude of the Z and W boson masses mz and rri\y, respectively, and also of the 
scalar masses; and the seesaw scale uir, much larger than mz, which is the scale of the 
singlet-neutrino masses. A small mass scale mo, below the electroweak scale, gives the 
order of magnitude of the elements of the Dirac neutrino mass matrix Mo, in such a way 
that the light-neutrino masses end up being of order m 2 D /m R . 

It was noticed in Ref. || that, in this framework, the dominant radiative corrections 
to the seesaw mechanism are quadratic in the Yukawa couplings, just like the tree-level 
masses, and that the radiative corrections are smaller than the tree-level results solely 
because of the appearance of the factor (167r 2 ) 1 in one-loop integrals. By using two 
Higgs doublets and only one right-handed neutrino singlet, this fact was later exploited 
in Ref. |§ for the radiative generation of the solar mass-squared difference; the ratio 
Amg/Am^ tm of solar over atmospheric mass-squared differences comes out of the correct 
order of magnitude for the large-mixing-angle MSW solution of the solar-neutrino prob- 
lem (for the main features of that solution, with Amg/Am^ tm being about 0.02, see for 
instance Ref. |4j]), with no other scale free to make adjustments. 

We want to extend here the calculations of Ref. ||, where the focus was the radiative 
generation of neutrino masses vanishing at tree level. In the present work we perform a 
general calculation of the dominant one-loop corrections to the seesaw mechanism, with 
emphasis on the demonstration of the finiteness and gauge-independence (in gauges) 
of those radiative corrections. As a result of this effort we shall present a formula, easily 
applicable in our extension of the SM, whose purpose is twofold: 

• To check the radiative stability of the seesaw mass matrix and of the relations 
derived therefrom, or else to compute the radiative corrections to those relations; 

• To explore the possibilities of radiatively generating Amg/Am^ tm . 
Unfortunately, the second point appears quite difficult to us, when one wants to combine 
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it with specific properties of the lepton mixing matrix U, like the small U e z and the 
practically maximal atmospheric mixing. 

This paper is organized as follows. In Section ^| we review the seesaw mechanism and 
our formalism for the multi-Higgs-doublet SM. We discuss the calculation of the neutrino 
self-energies and write down their one-loop integrals in Section [3]. The actual calculation 
of the dominant one-loop corrections to the seesaw mechanism and the demonstration 
of their finiteness and gauge-independence are carried out in Section |j. In Section |5| we 
give arguments in favour of the dominance of the corrections calculated in the previous 
section. In Section |6| we present the final results and a summary. 



2 Framework 

We consider the lepton sector of an extension of the Standard Model with rii left-handed 
doublets and rii right-handed charged singlets, i.e., til families, plus hr right-handed neu- 
trino fields and n# Higgs doublets. We denote the Higgs doublets by <pk {k — 1, 2, . . . , uh) 
and define <fik = iT2<j)%- The vacuum expectation value (VEV) of the neutral component 
of 4>k is Vf, I y/2. The Yukawa Lagrangian of the leptons is given by 

£y = - E {4^k + feA fc ) D L + H.c. (1) 
fc=i 

We employ a vector and matrix notation where £r, ur, and Dl are the vectors of the 
right-handed charged-lepton fields, of the right-handed neutrino singlets, and of the left- 
handed lepton doublets, respectively. The Yukawa coupling matrices T k are til X ul, 
while the are hr x til- The charged-lepton mass matrix Mg and the Dirac neutrino 
mass matrix are 

M e = -7nY. v l T k and M D = -^=J2 v ^k, (2) 
V2 k V2 k 

respectively. Without loss of generality, we assume Mi to be diagonal with real and 
positive diagonal elements: Mi = diag (m e , m^, m T , . . .). The mass terms for the neutrinos 
are 

- v R M D v L - - VrCMruI + H.c. , (3) 

where C is the charge-conjugation matrix and Mr is non-singular and symmetric. It is 
well known |10| that Eq. can be written in a compact form as a mass term with an 
(ni + tir) x (m + tir) symmetric mass matrix 




M D+M = „ " , (4) 



which may be diagonalized as 



U T M D+M U = rh = diag (m 1 ,m 2 ,..., m nL+nR ) , (5) 
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where the are real and non-negative. In order to implement the seesaw mechanism |J 
we assume that the elements of Mr> are of order win and those of Mr are of order rriR, 
with rap "C mjj. Then, the neutrino masses with i — 1, 2, . . . , Ul are of order m 2 D /mR, 
while those with i = nt + 1, • • • , til + nR are of order vtir. It is useful to decompose the 
(riL + tir) x {n L + ur) unitary matrix U as f|, [11 



U 



u L 

Ur 



(6) 



where the submatrix Ul is til X (til + and the submatrix Ur is x (n L +Ur). With 
these submatrices, the left- and right-handed neutrinos are written as linear superpositions 
of the Ul + ur physical Majorana neutrino fields xo 



vl = Ul'Jlx and v R = Ur'Jrx , 

where 7^ = (1 — 75) /2 and jr = (1 + 75) /2 are the projectors of chirality. 
The unitarity relations for U are fill 



(7) 



u L ul 



1 



n L 5 



u R u 



>R U R — Inn 5 u LVr — vn L xn R , 

From Eq. (0) we shall need the two relations 

U* L mUl = and U R mU ] L 
It follows from Eqs. (§) and (g) that 



0, 



and £/£E/l + LT^ZTr = 1 



n L +n R ■ 



M 



D ■ 



n L +n R 



UrM d = UrUrtuUI = (1 

The leptonic charged-current Lagrangian is 
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UlUl) MJ\ 



mU\ j . 



Cc 



y/2 11 



IlUlX + H.c. 



(9) 



(10) 



Til 



where g is the SU{2) gauge coupling constant. The interaction of the Z boson with the 
neutrinos is given by 



4c } = j- W [il (uiu L ) - lR (u T L ui 



4c 



(12) 



where c w is the cosine of the Weinberg angle. 

A full account of our formalism for the scalar sector of the multi-Higgs-doublet SM 
is given in Ref. | IT[| . We use the notation S® for the neutral-scalar mass eigenfields; 
the vectors b G <D™ H characterize each of those fields — for their precise definitions see 
Refs. 0, |IT |. The Yukawa couplings of the neutral scalars S® to the neutrinos are given 
by 



dp 



(s°) = -^T,s° b x [(u R A b u L + uIaJu r ) 1l + (uIa\Ur + ulAiui) lR 



X, 
(13) 
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with A b = J^k^k^k- The neutral Goldstone b is given by the vector bz 

with (b z ) k = w k /v || [TTJ . Here, v — (\ 
obtain 



V\\ + \V 2 \ H h |fn H l 

«0 



1/2 



2mw/g- We thus 

n (14) 

Similarly, the charged-scalar mass eigenfields are characterized by vectors a G C nff . 
The Yukawa couplings of the charged scalars are given by 



M 



£y { S± ) = E S aZ [iR (KUr) - 11 (T a U L )] x + H.c. 



(15) 



where A a = J2k a k^-k and T a = J2k a t^k- The charged Goldstone boson G = S^ w is 
given by the vector aw with (aw) k = Vk/v; therefore 
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M n and T 



m w 



Mi . 



m w 



(16) 



3 Neutrino self-energies 

The seesaw mechanism || tells us that, at tree level, the mass matrix of the light neutrinos 
is given by 

M^ cc = —MqM^Md , (17) 

where M D and Mr appear in the full mass matrix A4d+m of Eq. ([|). The one- loop 
corrections generate an ni xn^ submatrix 8 Ml in the upper left-hand corner of A4d+m, 
where at tree level there is a zero submatrix. The submatrices Mr, and Mr also receive 
one-loop corrections, denoted SM D and SM R , respectively. We shall then have, for the 
effective light-neutrino mass matrix at one loop, 

M v = MT C + 5 Ml - 5MIMr X M d - M^M^5M D + MIMr X 5M r Mr X M d . (18) 

One-loop corrections to the neutrino masses originate in the one-loop neutrino self- 
energy S(p), where p is the neutrino momentum. We use the decomposition 

S(p) = A L (p 2 ) Hl + A R (V) i lR + B L (p 2 ) 11 + B R (p 2 ) 1r . (19) 

The non-absorptive parts of Al }R are Hermitian, while those of B^r are the Hermitian 
conjugates of each other. Furthermore, since the neutrino field vector x consists of Ma- 
jorana fields, the self-energy must fulfil the consistency condition 

S(p) = C [£(-p)] T C- 1 , (20) 

and that consistency condition translates into 

A L = A T R , B L = Bl, and B R = B T R . (21) 

After computing the neutrino self-energies one must renormalize the neutrino fields: 

XL = (l + \ z^j xl , hence X R = U + \ A) Xr , (22) 
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since xr = CXl- The neutrino mass matrices must also be renormalized. Mass countert- 
erms are indicated by 5 C . The upper-left submatrix of A4d+m, which we call M L , has no 
counterterm since it vanishes at tree level: 5 c Mi = 0. The operation 5 C only acts on Mp 
and on Mr, with, from the second Eq. (g), 

5 C M D = -L E [(^) A * + v k (^ A *)] • (23) 
V2 k 

Then, 

= ° (24) 

The renormalized neutrino self-energy, in the basis where the tree-level neutrino mass 
matrix is diagonal, has the structure |12|, [13], |H | 

-iS r (p) = -iE(p) +i^(z L + 4) J*7l + * g + zT h)ilR 

5 c rh + - mz L + - zjra) 7 L - z (Vm + ~ + - z T L vii^ ^ R . (25) 

We indicate one- loop contributions to by (J 1-100 ?. There are tadpole diagrams indi- 
cated by £ tad P° le ; but we do not include them in S(p). 

Now we explain our strategy for calculating the dominant one-loop corrections to the 
seesaw mechanism: 

i. We calculate S(p) in the basis where the tree- level neutrino mass matrix is diagonal. 

ii. We are especially interested in B L (p 2 ), because 

SM L = 5 l - loop M L = U* L B L (0) U{ . (26) 
We have used the inverse of the relation (H) to transform back to the basis of M. d+m- 

iii. Denoting by A a typical coupling constant of the coupling matrices A&, the light- 
neutrino masses are of order A 2 — see Eq. ([H])- Therefore, we are allowed to evaluate 
B L at p 2 = in Eq. fl2§. Note that 5M L is identical with 5 1_loop M L , since 5 C M L = 0. 
For the same reason, 

£1-1oo PMl mugt be finita We ghall 

see this explicitly in the 
next section, where we shall also demonstrate the gauge invariance of S^^M^. 

iv. On the other hand, the one-loop corrections to Mo and to Mr are given by 

6 1 -***M D = U R B L {0)Ul and S^Mr = U R B L (0) V T R , (27) 

respectively. We shall argue in Section [5] that one-loop corrections to Mr are irrel- 
evant, while those to Mjj lead to subdominant corrections in M. v . 

First we consider the contribution of the one-loop graph with W to S(p): 



rrio 



x 



(27T) 

(ul) (u L ) £j rwiL + (u%) m)ti -fw-i^ , (28) 



6 



where iS^ v (k — p) is the propagator of a W ± with momentum k—p, and d is the dimension 



of space-time. Clearly, EL (p) only contributes to A L (p 2 ) and A R (p 2 ). 

Next we consider the contribution from the Z boson. This yields, in a general 
gauge parameterized by the gauge parameter £z, 



44 ■/ (2n) d ? F 



mi 



(uluiUulu.) 



{k-pf 



{k-pY 



(k — p) — m% friz \_{k — p) — ^z m \ (k — p) — m% 



.(29) 



Now we consider the contribution of the one- loop graph with S^- It is 



d d k 



£ 



(2n) a ^ (k - pf - ml^r k 2 - m 2 

4A a ),(r a ^. + (^rr)jAX) 



(30) 



When is the charged Goldstone boson G ± ) one uses Eqs. (|16D together with Eq. (|T0|) 
to derive 



—2 



2m^ 



1 



£ 



(2vr) d (k - p) 2 - ^m 2 w ^ k 2 - m\ 
mi (ul), d (U L ) £j + (uj:) ie (Ut) ej m,}. 

In similar fashion, the neutral scalars S® contribute (neglecting tadpoles) 



(31) 



1 



E 



x 



(27r) a r (fc-p) a -mg ^^-m? 
U R A b U L + UlA^U* R ). L (u R A b U L + U T L A T h U R 



(32) 



If the scalar S b is the neutral Goldstone boson G°, then, using Eqs. ([14]) and (|T0|), 

^ 2 ' A 1 



4m^ 



(2vr) d (fc-p) 2 -^m|?F- 
x (mC^C^ + UlUlm) u (mUlU L + U T L U\ 



mf 



rn 



(33) 



4 The computation of SMi 

In this section, we explicitly perform the calculation of SMl, defined in Eq. (p6|). 

It is clear from Eq. (|3~0"D and from the third Eq. @ that the exchange of 5' ± contributes 
neither to 5 1_loop M L nor to 5 1_1o ° p Mr; it only contributes to 5 l ' loop M D . We therefore have 

SM L = 5M L (Z) + 6M L (G°) + £ 5M L . (34) 

bytb Z 
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The Z boson part is given by 



5M L (Z) = i 



d d k 



44 J (2tt 



d 



+ 



k 2 



k 2 



k 2 — m 2 z m 2 z \k 2 — ^z m z k 2 — m 2 



The neutral Goldstone boson contributes 

,2 , d d k 



5M L [G 



(0 



^m 2 w J (27r) d k 2 - £ z m| L k 2 -rh 2 L ' 



L k 2 -rh 2 L ' 
(35) 

(36) 



Using 77% = c w m z and the relation 



m 



u* fc2m xj\ = u* _ 

k 2 — m 2 k 2 — m 2 



ul 



(37) 



which holds by virtue of U\fhU\ = 0, it is clear that the ^-dependent terms in Eqs. 
and (|36|) cancel and therefore 5Ml is gauge-invariant. One has 



5M L (Z) + 6M L (a ) 



d d k 



44 J (2-n f k 2 -m\ L 



1 m d - w 



m 



m 2 y I k 2 — m 



(3f 



The contribution to 5Ml from the exchange of all neutral scalars, except the neutral 
Goldstone boson, is 



5M L ($) 



d d k 



III 



2 J (2ir) d k 2 -m\ b R k 2 



in 



(39) 



Now we perform the actual calculation of the integrals. Let the dimension of space- 
time be d = 4 — 2e, with e — > 0. We define the divergent quantity 



k = -e~ l + 7 - ln(47r) - 1 
where 7 is Euler's constant. Then, we obtain 



SM L [S° b 
Similarly, 
SM L (Z,1) 



AlUtm [k + \nm z + 



327T 



2 u R 



In Tb 
r b - 1 



U R A b , with r b 



m 



m 



(40) 



(41) 



9 l 



1Qty 2 C: 



2 n 2 L 



Ut17i I k H h lnm 2 + 



In r z 

T Z - 1 



" 1 

m 



Ul, with r z = ^ T , (42) 



m 



is the contribution to 5 Ml (Z) + 8 Ml (G°) from the term with d in Eq. fl38|); the term 
with — m 2 /m| yields 



SM L (Z, 2) 



9' 



i'\V 



rz ~ 1 



(43) 
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In order to demonstrate the cancellation of infinities we need the orthogonality relation 
of the vectors b G C" H f§: 

E + (bz) j (b z ) k = . (44) 

b+b z 

Remembering that A& = J2k^k^k, we find that all terms independent of the boson 
masses — in particular the infinities — drop out in the sum J2b=ib z $Ml (S®) + 5Ml (Z,2). 
We may thus write 

5K ( S °) = 3^ A ^^T^ Ab ' (45) 

5M' L (Z,2) = -—^Utm 3 ^-Ul. (46) 
b4ir z my V rz — 1 

Also, because of the first Eq. (||), the terms in SMl (Z, 1) of Eq. (£42|) which are proportional 
to m (but not the term with mlnm 2 ), in particular the infinities, drop out. Thus, 

SM' L (Z, 1) = -4y Ulm Ul . (47) 

16-K z cf u rz — 1 

It is obvious that 5M' L (Z, 1) = -4<5Af£ (Z, 2) We thus we have the final result 

- E 3^ ^ ^ + ^ ^ . (48) 



Using Eq. (10), the last term of the previous relation reads 

-^^^7^. (49) 
b47i z m^ v rz — 1 

This expression, together with mw = gv/2, shows that the Z contribution to SMl is of 
the same order of magnitude as the neutral-scalar contributions. 



5 Dominance of 5 Mi among the radiative corrections 

We now want to give arguments for the dominance of 6 Ml in the radiative corrections to 
M.^ cc . The main point in our argumentation is the observation that the terms in 5Ml 
in Eq. (^) are smaller than the tree-level masses of the light neutrinos solely because of 
the factor (167T 2 ) -1 from the one-loop integrals. In this context, we have to remember 
that, in our framework, it is natural to assume all scalar masses to be of the order of the 
electroweak scale. 

The heavy-neutrino masses m nL+ i, . . . , m nL+nR (or the elements of Mr) are free pa- 
rameters of the theory; therefore, corrections to Mr are irrelevant for us. We stress that 
this contrasts to the light-neutrino masses mi, m,2, • • • , Tn nL , which are calculable, i.e., 
they are functions of the parameters of the theory; this is reflected in the finiteness and 
gauge-invariance of SMl- 

Concentrating on S 1 ' 1o ° p Ml> as given by Eq. fl2~T|), we note that both and S® con- 



tribute to it, see Eqs. (|30|) and (jjffi) , respectively. Checking Eq. fl29|) , we note the absence 
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of a Z contribution. In the integral of Eq. (p0[), for charged-scalar exchange, the large 
mass scale tur is absent; therefore, if Y is a typical Yukawa coupling, then we must have 
5 1-loop M D (5'„ ) ~ Y 2 me, where rri£ is a generic charged-lepton mass. With rri£ ~ mo, 
we have Y 2 m 2 D /mR as the order of magnitude of the charged-scalar contribution to M. v , 
which is suppressed as compared to the SMl contribution. Turning now to S® exchange, 
with Eqs. fl27D and (^) the result of the one-loop calculation is 

5 Hoop M fl (S° b ) = ^ A fe f/ L m f A; + In m 2 + 4 A b . (50) 

Again, we arrive at p Mjj{S^) ~ 7 2 mB. Therefore, the contributions to M. v of 
Eq. (TH) from 5 1_loop Mo are negligible when compared to the ones from 5 Ml- 



Some remarks concerning the wave function renormalization matrix zl are in order. 
Usually, that matrix is defined in such a way that the renormalized fields remain in the 
mass basis []12 , I3[ 14] . With Eqs. (fL7|) and fl26|) we do not adopt this convention: the mass 



matrix of the light neutrinos, A4 U , is not diagonal at order Y 2 . This is consistent with the 
fact that, at zeroth order in the Yukawa couplings, the light neutrinos are massless and, 
therefore, degenerate. This degeneracy is lifted only at order Y 2 , and a large rotation, 
with angles of order one, will be needed, in general, for the diagonalization of M. v . Since 
we assume that Zl is "small" as compared to one, we do not include the large rotation 
in this matrix. In the same vein, concerning the terms mzi + z\m in Eq. (JZ5]), we are 
allowed to set to zero the light-neutrino masses. Using the explicit form of zl for on-shell 
renormalization of the heavy neutrinos, as found, for instance, in Refs. |TT], and taking 
into account that the small elements of U are of order mo/mR, we find that the wave 
function renormalization terms of Eq. (|25"D contribute at order Y 2 m 2 D /mR or g 2 m 2 D /mR 
to Ai u , and then we do not need to consider them further. 

The tadpole graphs contribute only to 5Md- For their treatment and their relation to 
the 5 c Vk terms in Eq. (|23|), see the extensive discussion in Ref. |T5| . They have no impact 
on our discussion. 



6 Conclusions 

We have argued in the preceding section that the dominant radiative corrections to M.^ cc 
are given by Eq. (|26|), which has contributions only from neutral-scalar and Z exchange; 
equation (^8|) is the explicit form of Eq. (|26|) after carrying out the one-loop calculations. 
These radiative corrections are, like the tree- level mass term of Eq. ([H]), quadratic in the 
Yukawa couplings. In Eq. ((48|) — as compared to M.^ cc — only the terms of first order in 
I/trr are relevant (ttir is the scale of Mr). We may use the approximations 

U R ~{0,W), with W*M R W* ~m = diag(m nt+1> ...,m ni+ft J, (51) 

where W is a unitary ur x ur matrix whose elements are not suppressed by m^/mjj (m_o 
is the scale of Md). We thus have the final result 

M v = -MIMr X M d + 8 Ml , (52) 
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with 

^ z 32tt 2 \m ml) 

+ M l W * f i In 5^ W^Mp . (53) 

The sum over the neutral scalars includes only the physical ones. The second line in 
Eq. (|53]) is the Z contribution. Note that if accidentally M^ ee = -M^W* (1/m) W^M D 
vanishes, then, even for a single Higgs doublet, the one- loop Z and Higgs contributions 
to 8 Ml will in general be non-zero |TB . 



Having more than one Higgs doublet leads to flavour- changing neutral interactions. 
Neutral-scalar exchange at tree level then leads to processes like ^ — > e ± e + e~; however, 
such processes only restrict elements of the coupling matrices Tk, which do not occur in 
M. v of Eq. (|52l) . Moreover, with more than one Higgs doublet and at one-loop level, ra- 
diative decays like /i ± — > e 7 proceed via charged-scalar exchange. However, moderately 
large charged-scalar masses, still in the range of the electroweak scale, can render their 
decay rates compatible with experimental bounds |J. 

Working in the physical basis of the charged leptons, Eqs. (|52| ) and (^) determine 
the neutrino masses and mixings in the Standard Model with an arbitrary number of 
Higgs doublets. As we have shown, 5 Ml is finite and gauge- invariant. Moreover, we have 
argued that all other one-loop contributions to M. u are of order Y 2 m 2 D /mR (or g 2 m 2 D / 'uir) , 
where Y is a typical Yukawa coupling constant; consequently, 5Ml dominates the one- loop 
corrections. 

The fact that the order of magnitude of the one-loop corrections in Eq. ( |53"D is smaller 
than the tree-level mass term solely by a factor of order (167r 2 ) 1 In (mjj/mo), where m 
is a typical neutral-scalar mass, can be exploited for generating the ratio Am|,/Am 2 tm 
radiatively, such that the order of magnitude of that ratio is precisely the one adequate 
for the large- mixing-angle MSW solution of the solar- neutrino problem |J. However, the 
question of whether this framework also allows to explain the specific features of lepton 
mixing, namely the small element U e s of the mixing matrix, and an atmospheric mixing 
angle close to 45°, requires further investigation. 
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